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Abstract. We present five newly found galaxy clusters at z > 0.6 from the XMM Large-Scale Structure Survey (XMM- 
LSS). All five objects are extended X-ray sources in the XMM images. For three of them we have sufficient spectroscopically 
confirmed member galaxies that an estimate of the velocity dispersion is possible: XLSSC 001 at z = 0.613 and cr v = 867+^ 
km s"\ XLSSC 002 at z = 0.772 and cr v = 524+ 2 ^ km s^'and XLSSC 003 at z = 0.839 and a v = 780!^ 7 km s -1 . These 
three clusters have X-ray bolometric luminosities L x ~ 1 - 3 x 10 44 erg s _I and temperatures 2-4 keV, and consequently are 
O i 1 less massive than previously known clusters at similar redshifts, but nevertheless they follow the low redshift scaling relations 

1~ \ between L x , T and cr,., within the limits of the measurement errors. One of the clusters, XLSSC 004, is detected independently 

as an overdensity of galaxies of a colour R-z' = 1.4 that matches the redshift of the central galaxy z = 0.87, although it cannot 
unambiguously be confirmed by the spectroscopic observations alone. The highest redshift candidate cluster pertaining to this 
paper, XLSSC 005, is most likely a double cluster complex at a redshift around unity, associated with an extended X-ray source 
with probable substructure. 
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1. Introduction On the other hand, a much more efficient strategy is to 

perform cluster searches in X-rays: galaxy clusters are one of 

Statistical samples of distant galaxy clusters can serve as a the most power ful X-ray emitting objects, their X-ray emis- 

powerful tool to study the large-scale structure of the Universe. sion is of thermal origin (free . free emission of the electrons 

To obtain cluster samples deep enough and over a large area from the gas trapped int0 me deep gravitationa i potential well 

with optical or near- infrare d observations is qui te a difficult of the cluster an(t in geri eral, in hydrostatic equilibrium) and 

raskfAg. Gunn et al.[l986 Postman et al.[199 6| Gonzales et dep ends on the square of the gas density and so projection 

al.|200H Stanford et al[[997] Dickinson [1993. Projection ef- effects m not very important . Another advantage is the fact 

fects and the small field-of-view of the current NIR imagers that high galactic latimde fields of the x . ray sky are popu i ate d 

interplay to make relatively large contiguous surveys of distant mosdy by tWQ distinct kinds of object; AGNs or QSQs which 

clustets a dubious game. appear as point-like, and galaxy clusters, which are extended. 

It is straightforward to distinguish between the two classes us- 
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other hand, AGNs/QSOs are point-like sources and generally 
have harder power-law spectra. The spectral distinction, how- 
ever requires sufficient photon statistics and is prone to many 
assumptions and uncertainties. 

Ever since the first imaging X-ray observations, galaxy 
clusters were one of the main targets. This has led to flour- 
ishing X-ray based cluster surveys, initiated by the Einstein 
Medium Sensitivity Survey (Gioia et al. 1990 1 and then con- 
tinued with ROSAT both at low redshifts (e.g. Ebeling et al. 
[T998ll2000l De Grandi et al. lT999l Bohringer et al. 1998) and 
high redshifts (e.g. Rosati et al 1998, Jones et al. 1998 Romer 
et al. EUUUI Ebeling et aL EUOll Gioia et al. l2U0Tl Burke et al. 
120031 Mullis et al. l2003l for a review see Rosati et al. l200H . 
From these cluster surveys we have gained knowledge about 
individual objects, scaling relations, studies of the physics of 
the intra-cluster medium, the large-scale distribution, the clus- 
tering properties and the evolution of galaxy clusters. 

Nowadays cluster searches and surveys are brought to a 
new level with Chandra and XMM. Some of these projects are 
a natural continuation of existing surveys from the previous X- 
ray missions, following-up interesting objects and to get better 
insights on cluster physics, while other surveys use the publicly 
available archives of XMM and Chandra X-ray observations - 
this is really an inexpensive way to get statistically significant 
number of serendipitously detected clusters at all redshifts (e.g. 
Romer et al. l200T1 Schwope et al. 120031 Kim et al. l2004t . 

The increased sensitivity of XMM and Chandra X-ray ob- 
servatories, however, can also be exploited to survey large con- 
tiguous areas of the sky with relatively short exposure times. A 
fainter cluster/group population at intermediate redshifts that 
was not accessible to the previous X-ray missions will thus be 
detected. In return, this will give us a unique view of the cosmic 
web at the mass scale traced by this population. 

This is the approach of the XMM Large-Scale Structure 
Survey (XMM-LSS, Pierre et al.EOOlal 1200TB . The final ob- 
jective is to map a contiguous region of the sky and to study 
the large-scale distribution and the clustering properties of the 
matter traced by the galaxy clusters and the QSO/AGN popu- 
lation. XMM-LSS is based on XMM observations and a sub- 
sequent multi-wavelength follow-up programme. The target 
XMM-LSS survey geometry and depth were chosen such that 
to have a statistically significant number of clusters so that 
the two-point correlation function of clusters at < z < 0.5 
and 0.5 < z < 1 can be calculated with a precision reaching 
that of the ROS AT-ESO Flux Limited X-ray Survey (REFLEX, 
Collins et al. 2000 1 which covers a much larger area but only 
goes out to z < 0.2. 

To achieve this goal we need a good X-ray cluster detec- 
tion procedure. Using extensive XMM image simulations we 
have developed an efficient pipeline procedure that uses the full 
XMM instrument sensitivity and indeed allows us to detect ex- 
tended sources down to flux levels of 5 x 10~ 15 erg s _1 cm" 2 
for 10 ks exposures (Valtchanov et al. 2001 Refregier et al. 
2002 1. The successful detection is quite crucial because even 
the mere counting of clusters at different redshifts can be used 
to constrain cosmological parameters (e.g. Mathiesen & Evrard 
fl998l Henry 120001 Borgani et al. l200T1 Refregier et al. 2002 1. 
These cosmoslogical constraints that come from cluster abun- 



dance evolution are complementary and independent to those 
from cosmic microwave background and supernovae studies. 

Each detected candidate cluster from the pipeline, depend- 
ing on its estimated redshift by photometric redshifts, is pro- 
grammed to a spectroscopic follow-up (z < 1) or NIR obser- 
vations (z > 1). The spectroscopic follow-up of the first sam- 
ple of candidate clusters in the XMM-LSS at z < 1 was pro- 
grammed for observations on Las Campanas/Magellan and on 
the ESO/VLT telescopes. The subject of this paper is to present 
the first results for a sample of 5 clusters at z > 0.6 while 
the low redshift sample is presented elsewhere (Willis et al., 
in preparation). 

The paper is organised as follows: first we present the X-ray 
data reduction and source detection (Sect.|2j, then we describe 
in Sect.|3]the optical identification, spectroscopic target selec- 
tion procedure and the observations. In Sect.|4] we present the 
data analysis results from the spectroscopic and X-ray observa- 
tions. Next we discuss each individual object (Sect. |5} and we 
end up with the conclusions (Sect. |6j. Except where is men- 
tioned, we use ACDM cosmology (Hq = 70 km s _1 Mpc~', 
Q m = 0.3, Qa = 0.7) for all cosmologically dependent param- 
eters. All X-ray luminosities are bolometric. 



2. X-ray data 

The candidate clusters for the first spectroscopic run were cho- 
sen from all XMM-LSS pointings received by August 2002. 
This includes 15 AO-1 pointings of 10 ks exposure and an- 
other 15 exposures of 20 ks from the Guaranteed Time XMM 
Medium Deep Survey (XMDS). All observations were of good 
quality, except two fields that suffered from high background 
contamination affecting more than 50% of the exposure time. 

2.1. Data reduction 

A detailed description of the pipeline used in the XMM-LSS 
data reduction will be presented elsewhere (Pacaud et al., in 
preparation). Here we just briefly mention the main steps. The 
raw X-ray observations (ODFs) are reduced by the standard 
XMM Science Analysis System (XMM-SAS) tasks emchain 
and epchain for MOS and PN detectors respectively. High 
background periods, related to soft protons, are excluded from 
the event lists and raw photon images with a given pixel scale 
(2.5"/pixel in this case) in different energy bands are then cre- 
ated. Subsequently the raw images for each instrument are fil- 
tered using "a trous" (with holes) iterative wavelet technique 
with a Poisson noise model and a threshold of 10~ 4 (equiva- 
lent to 3.7<x in the Gaussian case) for the significant wavelet 
coefficients (Starck & Pierre IT9981 Starck et al. IT9981 . Each 
filtered image is then exposure corrected and a mask map that 
includes bad pixels, CCD gaps and non-exposed CCD regions 
(generally parts outside the field-of-view of the telescope) are 
created. Wavelet-filtered, exposure-corrected images for each 
instrument in a given energy band are added together to form a 
compound (MOS1+MOS2+PN) single band image to be used 
in the first stage of the detection procedure. 
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2.2. Source detection 

Clusters of galaxies are extended sources in X-ray images. 
Their detection and correct classification is not trivial because 
of various peculiarities of the X-ray observations: Poisson 
noise regime, varying PSF as a function of the off-axis angle 
and the energy, the vignetting effect, and the geometry of the 
detector (CCD gaps). We use images in the [0.5-2] keV en- 
ergy band which is well suited for clusters and groups (Scharf 
2002 1. The detection procedure is based on the prescription 
of Valtchanov et al. (120011 and has three stages: wavelet fil- 
tering (see the previous section), detection and measurements. 
The wavelet filtered image is fed to SExtractor (Bertin & 
Arnouts 1996) for detection and characterisation. The classifi- 
cation to extended (clusters) and point-like sources (AGNs or 
QSOs) is done by using three parameters: the half-light radius, 
the FWHM from a Gaussian fit to the source and SExtractor 
stellarity index adapted for the X-ray observations. The classifi- 
cation task is complicated by the fact that distant and faint clus- 
ters are not too different from a PSF. Moreover, at greater off- 
axis angles the PSF shape can be quite distorted, although the 
half-energy width does not change significantly in the [0.5 - 2] 
keV energy band. This shape distortion can lead to a wrong 
classification and that is why we have constrained the clus- 
ter detection up to off-axis distances not greater than 12'. This 
strategy gives very good results using simulated XMM images 
and can be used for selecting cluster candidates with an excel- 
lent success rate (Valtchanov et al. 2001 ). 

In total, from the first 3.5 deg 2 (30 XMM pointings) of 
XMM-LSS observations, we find some 55 X-ray cluster can- 
didates, which corresponds to ~ 15 clusters per sq.deg., in 
good agreement with the cosmological predictions (see e.g. 
Refregieret al. 120021 . 

3. Optical data 

3.1. Optical imaging 

Deep images from the CFH12k camera on the Canada-France- 
Hawaii Telescope (CFHT) in BVPJ from th e VIRMOS-VLT 
Deep Survey (VVDS: Le Fevre et al. l200T1 McCracken et al. 
2003 1 were used for optical identifications. The observations 
were processed by the Terapix team 1 to produce astrometricalry 
and photometrically calibrated images and to create object cat- 
alogs by means of SExtractor. The definition and magnitude 
limits of the VVDS may be found in the web pages of the 
VIRMOS consortium 2 . The X-ray cluster candidates were as- 
signed to NEAR (z < 0.5) and MID (0.5 < z < 1) samples by 
visual inspection of the optical images and also using photo- 
metric redshift estimates. We assign all the confident extended 
sources, that have no obvious optical counterpart in the VVDS 
images, in the DIST (z > 1) class for subsequent NIR follow- 
up. The visual scan of the data was also indispensable for re- 
moval of spurious extended sources introduced by instrumental 
effects. In addition to the CFH12k data, we obtained observa- 
tions at the CTIO 4m telescope in R and z' which allowed us to 

1 http://terapix.iap.fr 

2 http://www.astrsp-mrs.fr/virmos/vvds.htm 



search for a red sequence of galaxies at a given X-ray position 
(Andreon et al. l2004> . 

3.2. Spectroscopic target selection 

Our driving objective was to confirm the clusters and to mea- 
sure their redshifts. We chose the ESO-VLT/FORS2 instru- 
ment in MXU multi-object spectroscopy mode because of the 
liberty to place large number of slits with different sizes and 
orientations and so to optimise the target selection, especially 
in regions of high galaxy density. Using f ims (the FORS 
Instrumental Mask Simulator tool), we placed about 30 slits 
on average in each mask, possibly placing inclined and longer 
ones. With the remaining available slits we randomly sampled 
the X-ray point-like population presented in the same area or 
just filled the mask with objects from the field, when there were 
no colour selected or X-ray objects. 

We have developed a visualisation procedure to facil- 
itate the selection of spectroscopic targets, optimising the 
chances to pick up cluster member galaxies. It uses all the 
available information in the field, combining optical (one 
band or pseudo-colour VRI image) and X-ray images, colour- 
magnitude and colour-colour diagrams, and photometric red- 
shift peaks with the corresponding probabilities. In this "multi- 
parametric space", first we have taken the objects above the 
photometric limit, imposed by exposure time and observational 
strategy constraints. Secondly, we have looked for objects with 
a plausible colour (V-I or R-I) for ellipticals at the estimated 
photometric redshift and also used the pseudo-colour images 
for all candidates for which we have V, R and I observations. 
The utility of colour information is justified as, out to redshift 
of unity, the major population of cluster cores is made of ellip- 
ticals with an old stellar population (e.g. Dressier et al. 119971 
Postman et al. 1998 1. Moreover, cluster ellipticals at high red- 
shift show-up as an overdensity of red objects in VRI com- 
posite images, proving useful in the target selection process. 
Examples of RGB images for some of the candidate clusters 
can be found in the web pages of the XMM-LSS consortium 3 . 

3.3. Spectroscopic observations 

The ESO-VLT/FORS2 spectroscopic observations in MXU 
mode were performed on 9-12 October 2002. For the 3 allo- 
cated nights in "visitor" mode we selected some 12 candidates 
from the NEAR and MID samples. NEAR clusters were in- 
cluded in order to provide targets in case of a poor weather and 
will be presented elsewhere (Willis et al., in preparation). The 
observing log and object designations for the five MID clusters 
pertaining to this paper are shown in Table[T] 

We have used the holographic grism 600RI+19 with the 
GG435 blocking filter that gives a good response from about 
6000 to 8500 A and a dispersion of 1.64 A/pixel with the stan- 
dard instrument resolution. The spectral resolution for slits of 
1.4" width as measured from the arc or sky lines was « 6 
A FWHM. The high resolution of the grism (A/AA = 1000) 
makes subsequent sky subtraction easier. 

3 http://vela.astro.ulg.ac.be/themes/spatial/xmm/LSS 
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Table 1. ESO-VLT observing log. The exposure time is the 
number of masks times the time, and each mask is split into 
two exposures. 



ID" 


Cluster 


Seeing 


Exposure 






arc sec 


min. 




09-Oct-2002 






001 


XLSS J022457. 1-034853 


1.0-1.2 


2x60 


005 


XLSS J022709.7-041805 


0.7-1.1 


2x90 




10-Oct-2002 






004 


XLSS J022530.2-050713 


0.8-2.7 


90+60 




10-Oct-2002 






002 


XLSS J022532.5-035510 


0.7-1.2 


2x60 


003 


XLSS J022738.2-031757 


0.7-1.1 


2x60 



8 - In the following, all clusters are referred to via the reference 
XLSSC plus the identification number, e.g. XLSSC 001, etc. 

For each mask we split the total observing time into two ex- 
posures in order to remove cosmic ray contaminations. We have 
used the IRAF task imcombine with the crreject algorithm 
- appropriate for rejecting cosmic ray events even with two 
images. The spectral reduction was performed using the stan- 
dard tools from IRAF: zero level exposure, fiat-fielding, cos- 
mic ray removal, aperture selection and sky subtraction using 
the apextract package. For wavelength calibration, we have 
used He-Ne-Ar lamp exposures observed at the end of the cor- 
responding night. The wavelength calibration was performed 
with a Chebyshev polynomial of 3rd degree and the residuals 
were kept inside +1.0 A with an rms scatter ~ 0.3 - 0.5 A. A 
standard star was observed in the beginning and at the end of 
each night with the same instrumental configuration in order 
to remove the instrumental response and therefore to transform 
the spectra to relative flux units. 

4. Results 

4.1. Redshifts 

The sky subtracted, wavelength and flux calibrated spectra 
were used to derive the redshifts. First we obtained a redshift 
estimate based on the Call H+K doublet, when present, and 
then we used the cross-correlation technique as implemented 
in the RVSAO package (Kurtz & Mink 1998 1 using an ellipti- 
cal galaxy template spectrum (Kinney et al. I1996l >. To improve 
the cross-correlation signal, we have masked the wavelength 
ranges of the strongest sky lines, where residuals from the sky 
subtraction could occur, and also the region of a strong atmo- 
spheric absorption at 7750-7800 A. 

The redshift distributions in each field are shown in Fig.^ 
To illustrate the redshift space overdensities we have applied an 
adaptive kernel smoothing over the redshifts (e.g. Silverman 
1986| Pisani [1993 1. There is both a clear redshift grouping 
and a spatial grouping associated with the X-ray emission in 
XLSSC 001, XLSSC 002 and XLSSC 003. For these three 
cases only we show zoomed in an inset a redshift histogram 
of the overdensity with a fixed bin size. The derived mean red- 
shift and the velocity dispersion, by means of bi-weighted es- 
timators of location and scale (Beers et al. 11990b . are shown 



Table 2. Spectroscopic redshift measurements. N ga i shows the 
total number of galaxies used to derive the cluster redshift. The 
number in parenthesis shows the number of galaxies with [Oil] 
3727 A emission line. See the text for details on XLSSC 004 
and XLSSC 005 redshifts. 



Name 


Redshift 


km/s 




XLSSC 001 


n A19S+ 0006 
0.6128_ 0006 


867*8 


29(9) 


XLSSC 002* 


fl 7779+0.0001 
U. / / 000 3 


574+267 
JZ "^-116 


11(5) 


XLSSC 003 


fl OQS7 +0 0007 
U.8j8/_ .0007 


780^ 7 


17(5) 


XLSSC 004 


~ 0.87 




4(1) 


XLSSC 005 


~ 1.0 




7(5) 



(*) - taking only the objects having redshifts in [0.76, 0.78]. 

in Table |2j the quoted errors are lcr bias-accelerated bootstrap 
errors (Efron & Tibshirani ll986> . 

For XLSSC 004 and XLSSC 005, the peaks in the red- 
shift distribution over the whole FORS field do not corre- 
spond to any significant spatial clustering, and no peak can 
be unambiguously associated with the extended X-ray source. 
Nevertheless, we show in Table |2] the most plausible redshift 
based on measurements of galaxies spatially coincident with 
the cluster X-ray emission (see Figs.0and|8j. 

4.2. X-ray analysis 

Deriving detailed cluster physical characteristics from the X- 
ray observations is not the main objective of XMM-LSS. 
Indeed, with exposure times of the order of 10 ks only a small 
fraction of clusters are expected to have enough photon statis- 
tics to allow reliable measurements of the mean temperature, 
or to distinguish possible AGN contamination. Nevertheless, 
we have derived the temperature and luminosity for the clus- 
ters pertaining to this paper and, in some cases, these param- 
eters were well constrained. These results however must be 
taken with caution as in some cases possible AGN contami- 
nation cannot be ruled out. 

For each cluster an X-ray spectrum was extracted from a 
region large enough to include the cluster emission. A back- 
ground spectrum was taken from an adjacent annulus. We re- 
moved in advance the contribution of all other sources within 
the cluster and background regions. A photon redistribution 
matrix (RMF) and ancillary region file (ARF) were created 
using XMMSAS : rmfgen/arfgen, including corrections for bad 
pixels and detector geometry. Finally the spectra from the three 
instruments MOS1, MOS2 and PN were regrouped to have at 
least 25 counts per bin. The extracted spectra were used to de- 
rive the observed characteristics shown in Table [3] without any 
recourse to a reference model. 

To get the global cluster X-ray characteristics, the binned 
spectra in each instrument were fitted to a mekal model of 
thermal plasma emission with photo-electric absorption us- 
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0.2 0.4 0.6 0.8 1.0 1.2 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Redshift Redshift 




XLSSC 005 




0.2 0.4 0.6 0.8 1.0 

Redshift 



Fig. 1. Redshift distributions of the spectroscopically observed objects for all five clusters. The adaptive kernel density estimation 
(Silverman 1986 1 is shown as a continuous line over the redshift measurements. In cases where an overdensity in redshift space 
corresponds to a spatial correlation we show a histogram of the zoomed redshift region: for XLSSC 001 from 0.6 to 0.63, bin 
size 0.0027; XLSSC 002: 0.76 < z < 0.80, bin size 0.0025; XLSSC 003: 0.82 < z < 0.86, bin size 0.004. 



ing XSPEC (Arnaud fT996l also see the XSPEC manual 4 for the 
mekal model and corresponding references). The energy range 
used in the fit is [0.3 - 10] keV for MOS instruments, but for 
PN we used [0.3 - 7.5] keV in order to avoid an instrumental 



4 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/ 



emission feature at 8 - 9 keV. We have kept only two free pa- 
rameters in the fitting: the temperature and the normalisation. 
The mean Galactic absorption of Nh - 2.5 x 10 2l) cirT 2 (Dickey 
& Lockman 1990 1, the metal abundance of Z = O.3Z0 and 
the redshift were fixed. For parameter estimation, we used the 
Cash statistic, modified to allow background subtraction (see 
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Table 3. X-ray parameters. "Region" is the extraction radius in arcmin and in comoving Mpc in brackets. "/?5o" is the object's 
half-light radius in arcsec in the [0.5-2] keV band, note that the PSF half-light radius at 10' off-axis and at 1 keV is 9". The 
"Exposure" is the average of the weighted live time events in the extraction region. 



Name 


Redshift 


Region 




off-axis 


Counts 


Exposure 












[0.2-10] keV 


MOS1+MOS2+PN 






arcmin(Mpc) 


arcsec 


arcmin 


ALL(MOS 1 +MOS2+PN) 


ks 


XLSSC 001 


0.6128 


1.7(0.78) 


32.2 


5.8 


1812(472+467+873) 


2x13.3+7.9 


XLSSC 002 


0.7722 


1.1(0.60) 


16.6 


9.9 


633(177+152+304) 


2x13.3+7.9 


XLSSC 003 


0.8378 


1.1(0.65) 


23.6 


6.0 


575(136+128+311) 


2x11.5+8.3 


XLSSC 004 


0.87 


0.5(0.27) 




3.8 


233(40+83+110) 


2x20.8+16.6 


XLSSC 005 


1.0 


0.83(0.56) 


21.1 


8.5 


353(50+76+227) 


24.5+24.9+21.5 



the XSPEC manual). The results are shown in Table @] The er- 
rors on the temperature are lcr, the errors for the flux Fx and 
the luminosity Lx were calculated varying the normalisation 
parameter in its lcr confidence interval, while the temperature 
was kept fixed at its best fit value. The bolometric X-ray lumi- 
nosity Lx was calculated using the energy range [0.01-30] keV, 
constrained only by the instrument response matrices (gener- 
ally from ~ 0.05 to ~ 20 keV rest frame). 

For XLSSC 004 and XLSSC 005 photon statistics are insuf- 
ficient to allow any plausible parameter estimation so we have 
estimated the flux and the luminosity by fixing the temperature 
at 4 keV. The uncertainties in this case came from varying only 
the normalisation and thus can be quite underestimated. 

The cluster masses were estimated from the local L x - M 
relation (Reiprich & Bohringer 2002) assuming no evolution. 
The errors on the mass, reported in Table |4] only reflect the 
uncertainty on the luminosities and so are significantly under- 
estimated, given the large (up to an order of magnitude!) scatter 
in the local Lx - M relation. 

The X-ray spectra, mekal fit and the corresponding residu- 
als for XLSSC 001, XLSSC 002 and XLSSC 003 are shown 
in Fig. |3 and the derived model characteristics are given in 
Table|H 



Table 4. Derived, model dependent, cluster X-ray characteris- 
tics. The data presented for XLSSC 004 and XLSSC 005 is for 
fixed temperature of 4 keV and the uncertainties are probably 
underestimated. The reported masses should be treated with ex- 
treme caution as the scatter in the local Lx-M relation is large. 
See text for details. 



Name 


Redshift 


kT 
keV 


F a 

X 


Tj> 


M c 


XLSSC 001 


0.6128 


3 4 +L1 

J - -0.8 


8 9 +10 


3 5 +03 

J - -0.3 


A 4+0.2 


XLSSC 002 


0.7722 


2 2 +13 


2 6 +0 ' 7 


1 T 0A 

'-■ ' -0.3 


9 Q+0.5 


XLSSC 003 


0.8378 


4 r 47 


4 4+1-0 

^•^-0.8 


t 7+O.9 
■'• -0.6 




XLSSC 004 


0.87 


4 


0.23+°l 


21 +01 


0.8*» 


XLSSC 005 


1.0 


4 
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" - the unabsorbed flux F x in the [0.5-2] keV band and in units of 

10" 14 erg s" 1 cm" 2 . 

* - the luminosity L x is pseudo-bolometric ([0.01-30] keV) and in 
units of 10 44 erg s -1 , 

c - the total mass M in units of 10 14 Mq, based on the local L x - M 
relation. 



4.3. Scaling relations 

The scalings between the bolometric luminosity, temperature 
and velocity dispersion from the clusters in the compilation of 
Wu et al. ©99 1 are shown in Fig. [3] together with the ob- 
jects pertaining to this paper. XLSSC 001, XLSSC 002 and 
XLSSC 003 are in good agreement with the scaling relations 
L x -o- v ,L x -T and cr„ - T (Xue & Wu l^OTJUl Arnaud & Evrard 
I1999> . although caution must be taken because of possible bi- 
ases on the temperature introduced by the insufficient photon 
statistics (see e.g. Fairley et al. 2000 1. What is important to note 
is that the three clusters are at lower temperatures and luminosi- 
ties than the other known clusters at high redshift. In that sense 
we are starting to reach moderately massive systems at high 
redshift, whilst from ROSAT or EINSTEIN surveys the distant 
clusters are rather more luminous and consequently more mas- 
sive systems. 



5. Individual objects 

The optical/X-ray overlays for all candidates are shown in 
Figs l4l8l The optical images are from CFHT/CFH12k cam- 
era in the I-band, the X-ray contours are from wavelet fil- 
tered MOS 1 +MOS2+PN images (see Sec.0. For XLSSC 001 , 
XLSSC 002 and XLSSC 003, the objects from the inset his- 
tograms of their redshift distribution (Fig. [3 are denoted as 
boxes and triangles, the latter for the galaxies with [Oil] 3727 
A emission line. When there was no obvious redshift peak 
(XLSSC 004 and XLSSC 005) we indicate the redshift of each 
spectroscopically observed galaxy. 

Based on the morphology, redshift distribution and the vi- 
sual appearance, we define three broad cluster classes: relaxed, 
relaxing and collapsing. 
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Fig. 2. X-ray spectra, model fit and the corresponding residuals 
for MOS 1 (in black, the continuous histogram shows the model 
fit), MOS2 (in red, the dotted histogram is the model fit) and PN 
(in green, the model fit is the dashed histogram). The spectra 
were initially regrouped to have at least 25 photons per energy 
bin. 



5.1. Relaxed cluster: XLSSC 003 

XLSSC 003 (Fig. |6} has the appearance of a relaxed cluster 
with velocity distribution close to a Gaussian. Normality tests 
(Anderson-Darling and Shapiro-Wilk, see e.g. D'Agostino & 
Stephens 1986 1 accept the null hypothesis of the distribution 
being drawn from a normal distribution at 90% confidence, but 
this must be taken with caution as the number of objects is 
small. The X-ray emission is slightly elongated and there is no 
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Fig. 3. The scaling relations: luminosity-velocity dispersion 
(Lx - <t v , upper panel), luminosity-temperature (Lx - T, middle 
panel) and velocity dispersion-temperature (lower panel). The 
clusters from the compilation of Wu et al. ( 1999) are denoted 
as diamonds and those at z > 0.5 are shown as blue boxes with 
error bars. The best fit from Xue & Wu ( 2000 1 for Lx - cr v and 
cr,, - T and from Arnaud & Evrard ( 19991 for Lx-T are shown 
as lines. The objects pertaining to this paper are shown as filled 
red stars with the corresponding error bars. For this compari- 
son only, their luminosity was transformed to Einstein-de Sitter 
cosmology with Hq - 50 km s~' and A = (qo = 0.5). 
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evidence of substructure. The brightest cluster galaxy coincides 
with the X-ray emission peak and is at rest with respect to the 
cluster mean (zbcg - 0.8382, Av = 65 km s _1 ). 

As can be seen from Tabs.[3]and|3] the low photon statistics 
do not allow the temperature to be constrained, although if we 
assume an equipartition between the galaxy velocities and the 
gas, i.e. kT = (j,m p (r%, we derive a temperature of kT = 3.9 + 1.5 
keV - in excellent agreement with the best fit value. 

5.2. Relaxing cluster: XLSSC 001 

Morphologically XLSSC 001 (Fig. |4} looks like present day 
still relaxing clusters: there is an obvious centre formed by the 
brightest cluster galaxies displaced at ~ 10" from the X-ray 
emission peak. The brightest cluster galaxy (BCG) has a typ- 
ical giant elliptical galaxy absorption spectrum and is at 750 
km s _1 with respect to the mean cluster redshift. The velocity 
distribution is quite broad, but the null hypothesis that the ob- 
served distribution is drawn from normal cannot be ruled out 
- the Anderson-Darling statistics value is 0.527 which is less 
than the critical value 0.616 needed to reject the normality as- 
sumption at 90% confidence level. 

There are two other X-ray sources projected over the clus- 
ter emission. One is the subclump to the north-east, which is 
most likely associated with the cluster as there are three cluster 
members within the X-ray contours (two of which are emission 
line galaxies) and no obvious optical counterpart at the X-ray 
peak. Unfortunately, due to instrumental problems, we could 
not get the optical spectrum of the south-east point-like source. 

The X-ray characteristics of the cluster are relatively well 
constrained as we have about 1600 photons in the spectrum 
and the cluster is of low temperature. The low luminosity of 
the cluster is in agreement with its temperature and the velocity 
dispersion, as compared to the scaling relations of local clusters 
(Fig.0. 

5.3. Collapsing cluster: XLSSC 002 

XLSSC 002 (Fig. |5J resembles a point-like source in X-rays, 
although its half-light radius is R50 « 17" (see Table [3] to 
be compared to the PSF half-energy radius of 9") and the 
extension log-likelihood (as reported by the XMM-SAS task 
eraldetect) is 7.4, corresponding to probability of 6 x 10~ 4 of 
being a point-like source. The 8 centrally located galaxies form 
a crown-like figure, and five are in a very narrow redshift range 
0.771 - 0.775 most likely forming a compact core. The bright- 
est cluster galaxy is at redshift coinciding with the mean of the 
centrally located galaxies (zbcg = 0.7714) and has a typical 
giant elliptical galaxy absorption spectrum. 

We have to note that there are galaxies at similar redshifts 
(z » 0.77) all over the field as can be seen in Fig. [5] None of the 
other X-ray sources present in the field are cluster members. 

The temperature is quite tightly constrained by the X-ray 
spectral fit even though the photon statistics seems rather poor 
-480 photons in the [0.2-10] keV band in total. This is possible 
as the cluster is at low temperature. The relations between clus- 



ter global characteristics, Lx, T and o> (Fig. [3} are consistent 
with those observed from local clusters. 

There is a large fraction of emission-line galaxies in the 
centre (5 from 8) as identified by the presence of the [Oil] 
3727 A line. This line is not broad in Type 1 AGN, and thus 
it is difficult to assess the presence of an AGN with only this 
one emission feature - most of the diagnostic line ratios that 
can be used involve lines (H B , [OIII] 4959+5007 A, H a ) out- 
side the range of the grism RI600 used in this spectroscopic 
run. Optic al-to-X-ray flux ratios from deep surveys suggest 
that AGN with magnitudes typical of these galaxies have X-ray 
fluxes covering a wide range, between ^8xl0 -16 and ^8xl0 -15 
erg cirT 2 s _1 (0.5-2 keV), with the mean being «;2xl0" 15 erg 
cirT 2 s _1 (eg. Mainieri et al. 2002). Even in the unlikely event 
that all 5 emission line galaxies contain X-ray emitting AGN at 
this mean flux, the total would only be »40% of the measured 
cluster flux. It is perhaps more likely that most are Butcher- 
Oemler star-forming galaxies, implying an interestingly high 
star-forming fraction (eg. Fairley et al. 2002 find blue fractions 
up to 40% in low Lx clusters at z~0.3). Pending further studies 
of the galaxy population, we conclude that some contamina- 
tion by AGN may be possible, but that it is unlikely that all the 
X-ray flux originates in AGN. 

Based on the morphology and the redshift distribution, the 
most plausible interpretation is that the cluster is in its early 
stage of formation: a compact core is already formed and the 
accretion of matter from the nearby large-scale structure is un- 
derway - many galaxies in the field have redshifts around the 
redshift of the central core. 

5.4. Other cases 

XLSSC 004: The cluster X-ray emission (see Fig. lies be- 
tween two bright X-ray sources identified as QSOs at z — 1.12 
(the closest to the north-west, with 486 counts in the [0.5-5] 
keV band and V = 21.2, R = 20.9, I = 20.8) and at z = 1.19 
(to the south-east, with 378 counts and V = 18.5, R = 18.5, 1 
= 18.0). These two QSOs are located in diametrically opposed 
directions at 0.9' and 1 .2' from the centre of the cluster X-ray 
emission, respectively. Such a configuration is very unlikely 
to occur by chance. Our estimated surface density is approxi- 
mately 30 X-ray sources per deg 2 over a field of 2.8 deg 2 in the 
XMM-LSS survey with detected counts greater than 300, sub- 
sequently the probability of detecting two X-ray sources within 
2' from the cluster X-ray centre is about 1%. Given that these 
two quasars are located just behind a putative foreground clus- 
ter at z = 0.9 (see below), gravitational lensing could possibly 
affect the real positions and fluxes of these two background ob- 
jects. However, the hypothesis that these QSOs could be two 
images of a single source lensed by a foreground cluster does 
not hold as their optical spectra show a significant redshift dif- 
ference Az - 0.07. A more detailed discussion on this unusual 
association will be presented elsewhere (Jean et al., in prepara- 
tion). 

We cannot confirm the reality of the cluster from the spec- 
troscopic observations only - the brightest galaxy inside the X- 
ray emission is at z = 0.868 and there is another one at 15" to 
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the south-east at a similar redshift, z = 0.874. Two more galax- 
ies at north-west are at similar redshift 0.875, although outside 
the detected X-ray emission. 

XLSSC 004 is, however, independently detected on a scale 
of 1 arcmin diameter (marked as a thick green circle in Fig. 01 
as an overdensity of galaxies with a similar colour of R-z' = 
1.4 m , and cluster optical luminosity function compatible with 
z = 0.9 (Andreon et al. EOOH . 

The peak in the redshift distribution at z ~ 0.68 (see Fig.Q 
does not correspond to any spatial grouping. 

XLSSC 005: This is a very interesting and complicated 
case. The morphology of the X-ray emission suggests a bi- 
modality, with a central extended source and a possible sub- 
clump to the south-east (see Fig. [8} . It is not trivial, however, 
to disentangle the centre, where galaxies at z = 0.92 - 0.96 
and z = 1 .02 - 1 .06 are mixed. The redshift space distance is 
rather significant: Az ~ 0.1 (~ 15000 km s~') which rules out 
the galaxies belonging to one and the same cluster. 

This cluster has been observed in the K' band with ESO- 
NTT/SOFI in the framework of the XMM-LSS VIRMOS Deep 
Survey over an area of 0.25 deg 2 (Iovino et al. 2004, in prepa- 
ration). The exposure of 1 .5 hours resulted in a catalog down to 
a limiting magnitude K\g ~ 22.7. A search for overdensities 
in photometric redshift slices of Az p i lol = 0.2 was made, with 
Zphot calculated from BVRIK' photometry. As a result, there is 
a strong detection with signal-to-noise (S/N) ~ 5 of an overden- 
sity at z ~ 0.9 centred on the south-east extension of the X-ray 
emission. Observations in K at ESO-NTT/SOFI and R and z' 
at CTIO (Andreon et al. 2004) also indicate two overdensities 
of galaxies with colours compatible with z ~ 1 . 

The large fraction of emission line objects is also striking. 
From 5 objects near the X-ray peak with redshifts between 
0.92 - 1.06, 3 have detectable [OII]3727A emission in their 
spectra. As in the case of XLSSC 003, it is difficult to assess 
the presence of an AGN from this one emission line. At these 
redshifts, broad [MgII]2800A would be detectable, but in no 
case was it found. 

Based on the X-ray-to-optical flux ratios it is possible that 
some contribution to the overall X-ray emission might be due 
to AGN, but unlikely that all of it is due to AGN. 

Our preliminary interpretation is that we are observing two 
clusters in projection: one at z ~ 0.93, forming the main body 
of the extended X-ray source, and the more distant one at 
z = 1.02, that is partially covered. The redshifts of the spec- 
troscopically observed objects inside the X-ray emission span 
all of the interval from z ~ 0.92 to z ~ 1.05. This could be an 
indication of that some of these galaxies may lie on a filament 
connecting both clusters. Their number, however, is insufficient 
for any strong claim. 

Note that XLSSC 005 is very similar to the case of 
RXJ1053.7+5735 (Hasinger et al. [19^81 - a double cluster at 
z = 1.14 (Hashimoto et al. 2004), which was repeatedly re- 
ported in the literature at z — 1.26, mistakenly taken a cD- 
like galaxy from the background as the cluster central galaxy. 
Thanks to the deepest XMM single pointing observation it was 
possible to disentangle RXJ1053. 7+5735. Unfortunately we 
cannot apply the same technique as the exposure time in our 
case is ~ 50 times shorter. XLSSC 005 is however a good tar- 
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get for deep integral field spectroscopy observations (e.g.IFU 
of the VIMOS instrument) that will help us resolve this ex- 
tremely interesting case. 



Fig. 4. XLSSC 001. CFH12k I-band image of a 7' x 7' field is 
shown. North is up and east to the left. The X-ray contours run 
from 0. 1 to 5 counts/pixel (2.5" pixel size) with 10 levels in log 
space. Objects in the 0.6 < z < 0.63 redshift range are denoted 
by boxes and triangles (for emission-line objects). 



Fig. 5. XLSSC 002. AT xT field is shown and the objects 
at 0.76 < z < 0.81 are denoted by boxes and triangles (for 
emission-line objects). The X-ray contours run from 0.1 to 5 
photons/pixel with 10 levels in log space, with the lowest clus- 
ter contour at 0. 1 photons/pixel. 



Fig. 6. XLSSC 003. Only the central 5' x5' region of the field is 
shown and the objects in 0.82 < z < 0.86 are denoted by boxes 
and triangles (for emission-line objects). The X-ray contours 
run from 0. 1 to 5 photons/pixel with 10 levels in log space, and 
the lowest cluster contour is at 0. 1 photons/pixel. 



Fig. 7. XLSSC 004. Only the central 3.4' x 3.4' region of the 
field is shown and the redshifts for the spectroscopically ob- 
served objects at z > 0.8 are indicated. The triangles denote 
emission-line objects. The X-ray contours run from 0.15 to 5 
photons/pixel with 10 levels in log space, with the lowest clus- 
ter contour at 0.15 photons/pixel. The optical detection from 
CTIO Rz' data is shown as a green circle. 



Fig. 8. XLSSC 005. Only the central 3' x 3' region of the 
field is shown and the redshifts for the spectroscopically ob- 
served objects at z > 0.2 are indicated. The triangles denote 
emission-line objects. The X-ray contours run from 0.15 to 5 
photons/pixel with 10 levels in log space, with the lowest clus- 
ter contour at 0.15 photons/pixel. 



6. Conclusions 

We present five newly discovered high redshift, X-ray selected 
clusters in the XMM Large-Scale Structure Survey. The detec- 
tion and classification in X-rays, the subsequent optical identi- 
fication and spectroscopic target selection demonstrate the effi- 
ciency of the programme and the associated follow-up. 
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6.1. Highlights 

- The five newly discovered X-ray clusters at z > 0.6 were 
observed at ESO-VLT for 11.5 hours in total. For three of 
them, we have measured more than 15 galaxies with con- 
cordant redshifts and thus obtained a viable estimate of the 
velocity dispersion. 

- Thanks to the optimised target selection, we were able to 
get an estimate of the velocity dispersion of XLSSC 003 
at z = 0.84 in a reasonable VLT time of 2 hours. To show 
the potential of the XMM-LSS it is interesting to note that 
to date there are 5 known clusters in the literature (two of 
them are X-ray selected) at z > 0.8 with more than 10 spec- 
troscopically measured concordant redshifts. 

- We have detected a complex structure at a redshift of unity 
- XLSSC 005. The X-ray morphology and the presence of 
galaxy overdensity in photometric redshift space and in K, 
R and z' colour space allow us to speculate that most likely 
we see two clusters in projection. Some galaxies included 
in the spectroscopic observations may reside in a possible 
filament connecting the two clusters. 

- From CTIO Rz' observations, XLSSC 004 is detected as an 
overdensity of R-z' = 1.4'" galaxies, matching the redshift 
z ~ 0.9 of two bright and centrally located galaxies inside 
the X-ray extended emission. 

- From Tables [3] and @] and Fig. [3] it is obvious that for the 
first time we are unveiling, in a systematic manner, moder- 
ate mass clusters in the redshift range 0.5 < z < 1. Most of 
these objects are weak, extended X-ray sources and conse- 
quently would have been difficult to classify in typical deep 
ROSAT/PSPC pointings because of the insufficient pho- 
ton statistics in addition to the worse PSF, although there 
are successful ROSAT identifications of a handful of clus- 
ters at similar redshifts and mass scales (eg. Maughan et 
al. l2003l Vikhlinin et al. l2002l Thus, with XMM-LSS we 
are starting to fill the cluster database with significant num- 
ber of objects at high redshift from the middle of the mass 
function. This is a great improvement upon ROSAT based 
cluster surveys and will bring important new information 
on cosmological constraints as well as on non-gravitational 
physics in clusters (as the effects like pre-heating and feed- 
back are presumably more important in lower mass sys- 
tems). 

6.2. Prospects 

In the near future, our goal will be to confirm, via spectroscopic 
identifications, all XMM-LSS clusters down to ~ 8 x 10~ 15 erg 
s _1 cirT 2 over ~ 8 sq.deg. of the AO-1 and AO-2 observations. 
This will form a complete sample of about 100 X-ray selected 
clusters at < z < 1 . A follow-up of the most interesting ob- 
jects at z ~ 1 is also previewed with integral-field spectroscopy 
and possibly deeper targeted X-ray observations with XMM or 
Chandra for better insight into the cluster physics. 
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